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The  current  meter  records  collected  at  three  sites  in 
the  Gulf  of  Mexico  during  the  passage  of  Hurricane  Frederic 
are  analyzed  to  determine  the  storm-induced  flow  at  various 
ocean  depths,  determine  the  associated  =aergy  increase  and 
decay,  and  compare  these  observations  to  similar  results 
from  a  numerical  model.  The  records  at  the  two  deeper  sites 
are  rather  unigue  because  they  are  within  100  km  of  the  hur- 
ricane track.  Pre-storm  conditions  are  controlled  by  topog- 
raphy, and  as  the  storm  passes  there  is  an  abrupt  change  in 
the  direction  of  flow  and  initiation  of  a  strong  inertial 
response  at  all  levels  of  the  two  deeper  sites.  After  this 
initial  surge,  the  residual  flow  tends  toward  the  pre-storm 
direction.  The  horizontal  kinetic  energy  associated  with 
inertial  motion  is  calculated.  The  energy  increase  and 
decay   is    shown   to   vary   with  depth. 

An  embedded  mixed-layer  ocean  circulation  model  (Adamec 
et  al,  1981)  is  forced  with  an  idealized  storm  translating 
at  the  same  speed  (7.5m  s-1)  as  Fredsric.  The  abrupt 
response  and  strong  inertial  component  predicted  by  the 
model   is    qualitatively    similar   to   the    obs=rvations. 
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I-       III  REDUCTION 

Several  studies  of  the  response  of  the  ocean  to  hurri- 
cane passage  have  been  made  by  such  authors  as  Fisher 
(1958),  Leipper  (1967),  Wright  (1969),  Pudov  et  al  (1979) 
and   Fedorov      et    al    (1979).  A    review   of      the   hydrographic, 

i.e.  sea-surface  temperature  (SST),  salinity  temperature 
depth  (STD),  and  expendable  bathythermograph  (XBT)  ,  sur- 
veys of  several  of  these  authsrs  is  given  in  Table  I  from 
Price    (1981).  The    majority   of    these      observations   concen- 

trate on  the  asymmetrical  SST  response,  and  the  expected 
ocean  mixed-layer  response  to  hurricane  passage  of  vertical 
mixing  as  the  hurricane  approaches,  followed  by  upwelling  in 
the    wake    of      the   storm    (see  e.c[.  Leipper,       1967;      Friese, 

1977).  Additionally,  oscillations  on  the  order  of  the  local 
inertial  period  have  been  noted  in  the  temperature  and 
current  fields  in  the  wake  of  open-ocean  storms.  The  iner- 
tial response  becomes  complicated  in  the  shoal  waters  of 
continental  shelves  and  the  ability  to  detect  inertial 
motion  is  severely  restricted  (flayer  et  al,  1981).  The  hur- 
ricane-induced upwelling  decreases  the  nixed  layer  depth. 
The  depth  reaches  a  local  minimum  about  one  half  inertial 
period  after   eye   passage  and   continues   to    oscillate   at    about 
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the  local  inertial  period  for  some  period  after  the  storm 
passage.  Geisler  (1970)  presented  the  relationship  between 
hurricane  translation  speed  and  the  oscillatory  response  of 
the  thermocline.  If  the  hurricane  translation  speed  is 
greater  than  the  internal  phase  speed,  typically  about  2  m 
rlt  and  if  its  horizontal  scale  is  comparable  to  the  inter- 
nal Eossby  radius,  currents  throughout  the  affected  area  are 
controlled  by  a  balance  between  the  centrifugal  and  Coriolis 
accelerations. 

TABLE  I 

Hydrographic  Studies  of. the  Sea  Surface  Temperature  Response 

to  Hurricanes  (Price,  1931): 


Study:  Hurricane 


Method:  Region 


Hurricane 
Average       central 

UH  pressure      A  SST„ 

(in  s"')  (mb)  ("O 


Position  of*  &  SST- 


Leipper  ( 1967):  Hilda  (1964) 

extensive  post-hurricane 

hydrographic  survey:  Gulf  of 
Mexico 

3 

930 

-6 

pattern  is  generally 
unclear,  may  be  50  km 
to  left  of  track  (Fig.  8) 

Fedorov  «  a/.  (1979): 

Extensive  pre-  and  post- hurricane 

6 

980 

-2 

30  km  to  right  (Fig.  3, 

Ella  (1968) 

XBT  survey:  mid-Atlantic 

same  as  this  Fig.  la) 

Pudov  el  at.  (1979):  Tess 
(1975) 

extensive  post- hurricane  STD  survey: 
mid-Pacific 

6 

940 

-4 

75  km  to  right  (Fig.  1. 
same  as  this  Fig.  2a) 

Wright  ( 1969):  Shirley  ( 1965) 

1  pre-,  1  post-hurricane  XBT  section: 
vicinity  of  the  Kuroshio 

13 

935 

-3 

20  km  to  right  (Fig.  4) 

Jordan  (1964): 
Wanda  (1956) 
Clara  (1955) 

extensive  pre-  and  post-hurricane 
SST  reports  from  ships  of 
opportunity:  mid-Pacific 

16 

18 

920 
915 

-2 
-1 

150  km  to  right 
50  km  to  right 
(Figs.  2  and  3) 

^Estimates   made  by    Price    (19  81)     from   their   figures    noted 


The   response   to   Hurricane    Eloise      as    her    eye    passed   over 
the    National      Oceanic  and    Atmospheric      Administration's    Data 
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Buoy  Office  (NDBO)  EB-10  buoy  has  been  studied  extensively 
(see  e.j.  Martin,  1982;  Price,  1981;  Black  and  Withee, 
1976).  EB-10  was  located  in  the  central  Gulf  of  Mexico  and 
collected  the  first  open-ocean  data  under  a  hurricane  as 
reported  by  Withee  and  Johnson  (1976)  .  Mayer  et  al  (1981) 
reported  a  study  of  the  passage  of  Hurricane  Belle  over  the 
continental  shelf  of  the  New  York  Bight  on  10  August  1976. 
These  continental  shelf  responses  differ  from  the  open-ocean 
responses  due  to  the  topographic  influences  of  the  shelf  and 
large  gradients  in  the  physical  properties  between  the  shelf 
water  and   the  deep   ocean   water. 

Hurricane      Frederic   passed      through   the      Gulf   of      Mexico 
between    11      and    13      September,       1979.  The    satellite      data 

depicted  in  Figs.  1,  2,  and  3  together  with  whatever  XBT 
data  collected  during  hurricane  passage  are  usually  the  only 
type  of  data  available  to  study  these  geophysically, 
socially  and  economically  important  events.  However,  during 
this  period  the  Naval  Oceanographic  Office  (NAVOCEANO) ,  NSTL 
Station,  Bay  St.  Louis,  Mississippi  operated  three  data 
buoys  which  collected  a  unique  set  of  ocean  temperature  and 
current      data.  This      thesis      studies      the      ocean      current 

response    to  the      passage   of  Hurricane   Frederic      as   it    passed 
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near  ths  three  data  buoys.  Additionally  a  comparison  of  the 
real  data  to  the  results  provided  by  ths  three-dimensional 
ocean  model  of  Adamec  et  al  (198  1)  is  made.  The  hurricane 
forcing  of  the  model  is  idealized  rathsr  than  being  formu- 
lated to  represent  the  actual  hurricane,  but  a  translation 
speed  of    7.5    m   s~l,    the    same   as    for    Frederic,    is   used. 

The  buoy  current  data  provided  by  NAVOCEANO  indicate 
that  the  flow  associated  with  the  passage  of  Hurricane 
Frederic  had  a  large  inertial  component  after  the  storm 
passed.  The  rate  at  which  the  inertial  flow  damps,  and  the 
rate  at  which  it  propagates  with  depth,  are  determined  from 
the    tuoy    data   and   compared   to   the   model. 
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II.   DATA  PROCEDURES  AND  RESULTS 

A.   EACKGROUND 

Ten  current  meter  array  records  collected  at  three  sites 
on  the  continental  shelf  of  the  southeastern  United  States 
(Fig.  4)  during  the  passage  of  Hurricaie  Frederic  were  pro- 
vided by  NAVOCEANO.  As  depicted  in  Fig.  5,  the  arrays  con- 
sisted of  Aanderaa  spar  buoys  anchored  tD  the  bottom.  The 
anchor  cable  held  instruments  at  three  depths  at  each  of  the 
stations  one  and  two  and  at  four  depths  at  station  three. 
The  shallowest  instrument  depth  was  19  m  at  station  two, 
while  the  deepest  was  457  m  at  station  three.  The  bathyme- 
try (Fig.  5)  shows  all  three  stations  are  or.  the  continental 
shelf.  Data  were  provided  for  a  continuous  period  from  a 
few  days  before  the  storm  passage,  about  2100  Greenwich  Mean 
Time  (GMT)  12  September  1979,  until  a  few  days  after  the 
storm  passage. 

tfater  depth  at  station  one  is  only  about.  100  meters. 
Price  (1981)  reported  difficulties  with  aeters  in  a  similar 
coastal  environment  because  of  the  presence  of  very  strong 
horizontal  gradients   and  strong  topographic  effects   on  the 
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Figure  u.      3est    position   track   of   Hurricane    Frederic  and 
isitions   of   current    in 
K.     Shay,     NAV0CEAN3)   . 


positions    of^current^  meter  arrays    (provided    by 


currents.  Mayer  et  al  reported  intense,  first-mode,  near- 
inertial  period  oscillations  in  water  of  about  70  meter 
depth.  Similar  observations  at  sites  in  shallower  water  of 
about  50  meter  depth  resulted  in  only  weak,  heavily  damped, 
second-mode   oscillations. 

The    current    meter    records    used      in   this    thesis   also   show 
a    variation      in    response   among      the   stations      with    different 
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bottom  depths.  At  station  three  where  the  water  depth  is 
about  465  meters,  an  increase  in  current  nagnitude  starting 
about  1920  GMT  12  September  1979  is  very  apparent  in  the  raw 
data      record         (Fig-    6a).  This        record      also        shows      an 

oscillation  which  cycles  through  one  wavelength  about  once 
every     24    hours      (1440    minutes).  This   oscillation      decays 

smoothly  over  about  a  week,  and  at  first  glance  it  appears 
to  have  a  period  of  about  22  hours.  A  more  detailed  analy- 
sis cf  the  record  shows  the  zero  down-orossings  to  be  very 
nearly  1450  minutes  apart  during  the  largest  oscillation  in 
the  record.  The  average  period  for  seven  cycles  during  this 
interval  is  1430  minutes  (23.3  hours)  between  down-cross- 
ings- The  local  inertial  period  at  station  three  is  1473 
minutes-  Shortly  after  the  surface  response,  the  rapid 
increase  in  current  magnitude  is  also  obvious  at  the  depth 
of  457  meters  at  station  three  (see  Fig.  6b).  The  average 
period  of  the  oscillation  at  this  greater  depth,  1370  min- 
utes or  22-8  hours,  is  less  than  the  period  near  the  sur- 
face, but  the  damping  of  the  oscillation  is  not  obvious  at 
457  meters.  At  station  one  where  the  bottom  depth  is  only 
about  100  meters,  the  increase  in  current  speed  is  detecta- 
ble   in   the   raw   data    record      near    the    surface     (Fig.    7a)  ,      bat 
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any  period  of  oscillation  on  the  order  of  24  hours  or  a 
damping  rate  is  difficult  to  detect.  An  increase  in  current 
velocity  on  12  September  1979  near  the  bottom  at  station  one 
(Fig.  7b)  is  not  apparent  and  it  can  be  seen  that  current 
magnitudes  are  larger  on  9  September  than  they  are  on  12 
September. 

To  determine  if  the  observed  currents  were  extraordi- 
nary, some  idea  of  the  likely  currents  in  the  area  of  obser- 
vations is  helpful.  A  dominant  feature  of  the  Sulf  of 
Mexico  surface  currents  is  the  Loop  Current  (Fig.  3).  This 
is  a  clockwise  current  of  about  50  to  203  cm  s_l  and  90  to 
150  km  width  (Leipper,  1970)  which  enters  the  Gulf  in  the 
west  as  the  Yucatan  Current,  and  exits  through  the  Florida 
Straits  as  part  of  the  Gulf  Stream  system.  The  Yucatan 
Current  flows  north  from  Honduras  between  the  Yucatan 
Peninsula  and  Cuba  into  the  central  east  Gulf  and  forms  the 
western  section  of  the  Loop.  Further  north  the  current 
flows  east,  southeast,  and  east  again.  Elsewhere  in  the 
northeastern  Gulf,  the  surface  flow  is  generally  cyclonic 
along  the  coast  until  turning  southward  near  Louisiana. 
Eddies  are  known  to  have  detached  from  the  Loop  Current 
(Leipper,    1970)    and  could   move   into   the  area   of   the 
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Figure   6.      (J    components    of    the   raw   currant    meter   records    fo: 
7-21    September    at    iaoths   of    (a)     21    m    and    (b)     u57 
m    at    station    3. 
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Figure  7.   7  components  of  the  raw  currant  meter  record  for 
7-21  September  at  depths  of  (a)  49  m  and  (b)  92  m 

at  station  1. 
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instrument  arrays-  Leipper  (197  0)  shows  a  systematic  devel- 
opment and  breakdown  of  the  Loop  Current  which  places  the 
northern  edge  of  the  current  further  north  in  the  spring 
than  in  other  seasons  (Fig.  9),  Molinari  (1978)  suggests 
the  northward  intrusion  of  the  Loop  Curreat  is  not  seasonal, 
but  that  climatological  results  are  biased  by  temporal  sam- 
pling techniques,  Hurlburt  and  Thompson  (1980)  presented  a 
numerical  study  of  Loop  Current  intrusions  and  eddy  shed- 
ding. They  obtained  theoretical  expressions  for  the  eddy 
diameter  and  penetration  distance  of  the  Loop  Current  into 
the  Gulf  and  concluded  that  if  only  vorticity  dynamics  are 
considered,  the  interaction  between  the  Florida  Shelf  topog- 
raphy and  the  pressure  field  results  in  a  balance  which 
stops  the  northward  penetration  of  the  Loop  Current.  All  of 
the  current  array  stations  are  further  north  than  29  deg 
North,  and  although  the  waters  of  the  Loop  rarely  extend 
this  far  north,  Huh  et  al  (198  1)  examined  an  intrusion  of 
Loop  waters  as  far  north  as  the  data  buoys  which  collected 
the  raw  data  used  in  this  thesis.  It  is  not  possible  to 
tell  a  priori  whether  or  not  the  Loop  current  is  in  the  area 
of  the  current  meter  arrays  during  the  period  of 
observations,  but  it  seems  unlikely. 
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Figure  8.  Typical  surface  currents  in  the  Gulf  of  Mexico 
for  the  months  of  July,  August,  September.  Th 
general  location  of  the  buoy  arrays  is  indicat 
by   the   box.     (U.S.    Navy,    1965). 


Another  factor  which  might  be  expected  is  the  influence 
of  topography-  The  principle  cf  conservation  of  potential 
vorticity  causes  the  current  to  follow  the  ba-hymetric  con- 
tours. The  proximity  cf  the  buoys  to  the  coast  suggests 
that  near-shore  property  gradients  between  the  shelf  water 
and  the  deep  water  may  have  important  dynamic  and  thermody- 
namic consequences.  The  DeScto  Canyon  is  also  in  the  area 
cf  observations  and  may  tend  to  enhance  cross-shelf 
circ  ulation. 
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Figure  9.      The   spring   intrusion    of   1966    as    indicated   by 

overlays   of  the    150-meter   contour   lines   from   the 
22°C   topographies   of    all  spring    cruises   of    1966. 
Say    1964    and    November    and    December    1965   are 
included   for    supplemental    information    (Leioper, 
1970) . 

B.        RAW    DATA 

The      NAVOCEANO    data      consisted   of      north-south    (v)         and 

east-west    (u)    velocity   components  and  temperature    (T)       at    10 

minute      intervals.        Temperatures     above      21.5      °C   were      not 

recorded    due    to    thermistor    limitations.        Temperatures    above 
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21.5  °C  occurred  only  near  the  surface  and  complete  tempera- 
ture records  were  available  only  from  the  deeper 
instruments.  Also,  NAV0CEAN3  reported  apparent  internal 
clock  synchronization  problems  in  the  current  meter  record 
at   64   meters   depth,    station   on*. 

The  hurricane  passed  the  data  buoys  with  a  translation 
speed  of  about  7.5  m  s-1,  and  as  depicted  in  Fig.  4,  the  eye 
passed  about  100  km  to  the  west  of  station  three  at  2200  GMT 
12  September  1979  (Julian  day  255).  Maximum  winds  of  115 
Jets  occurred  at  a  radius  of  about  30  km  from  the  storm  cen- 
ter. The  inertial  period's  for  the  data  buoy  locations  range 
from  23  hours,  57  minutes  at  station  one  to  24  hours,  33 
minutes  at  station  three.  That  is,  the  inertial  period  at 
all    of      the  stations   is      approximately   one    day.  Since  the 

diurnal  tidal  period  is  24.8  hours  and  the  semi-diurnal 
tidal  period  is  12.4  hours,  difficulties  could  arise  in  sep- 
arating any  inertial  motion  from  tidal  motion  before  the 
storm  arrival.  A  spectral  analysis  would  most  likely  not 
have  sufficient  resolution  to  separate  inertial  and  diurnal 
motion.  The  first  harmonic  of  the  inertial  frequency  would 
also  be  indistinguishable  from  the  semi-iiurnai  tidal  fre- 
quency.        However,    the   large   increase   in    currents,      and   thus 
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kinetic  energy,   observable  in  the  raw  data  records  would  be 
easily  detected  in  an  energy  spectrum. 

The  initial  data  records  provided  by  NAVOCEANO  span 
non-coincident  periods.  However,  all  records  contained  data 
from  1920  GMT  07  September  1979  to  1820  GMT  21  September 
1979.  The  raw  data  records  (Figs.  6  and  7)  show  an  obvious 
response  to  the  hurricane  forcing  starting  about  12 
September  1979.  This  hurricane  forcing  results  in  increased 
currents,  with  maximum  values  of  130  cm  s~l  near  the  sur- 
face. Comparison  of  Figs.  6a  and  5b  shows  that  the  forcing 
response  was  transmitted  as  deep  as  457  neters  in  less  than 
one  half  an  inertial  period.  The  response  to  hurricane  pas- 
sage is  more  difficult  to  detect  in  most  of  the  raw  data 
records  of  stations  one  and  two  and  it  is  therefore  not  pos- 
sible to  estimate  the  rate  at  which  the  response  is  verti- 
cally propagated  with  depth  merely  by  looking  at  the  raw 
data  record.  The  ability  to  detect  the  increase  of  energy 
associated  with  storm  passage  at  any  station  and  at  any 
level  is  discussed  in  the  section  titled  "Procedures  for 
determining  relative  energies".  Although  this  thesis  does 
not  analyze  the  temperature  records  directly,  the  response 
to  hurricane  passage  is   also  apparent   in  the   temperature 
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record      (Fig.    10).  The      temperature        record      shows        an 

oscillation  in  the  thermocline  and  an  apparent  cooling  at 
the  437  in  depth  at  station  three,  followed  by  a  gradual 
warming   to   the   end   of  the   period. 
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Figure   10.      Raw    temperature   record   for   7-2  1    September    1979 
at  a   depth   of    437   on    at   station    3. 


C.        EATA    MANIPULATION 

The  initial  records  were  shortsned  to  include  only  the 
period  1920  GMT  7  September  1979  to  1820  GMT  18  September 
1979.  This  period  was  the  longest  possible  period  included 
by  all  records  and  was  a  logics. 1  choice  since  the  hurricane 
eye    passed   all  three    buoys    between   2100    3J1T    12    September    and 
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0300  GMT  13  September  1979.  By  examining  plots  of  the  data 
(e.3.  Figs.  6  and  7)  it  was  determined  that  significant 
storm  effects  with  some  pre-storm  and  soma  post-storm  data 
would   be    included   in   the   shortened  record. 

Progressive  vector  diagrams  (PVD's)  were  plotted  for 
each  record  during  the  above  period.  These  diagrams  indi- 
cated that  the  inertial  motion  was  superposed  on  a  mean 
current  which  was  quite  different  at  each  location.  The 
mean  current  was  not  of  primary  interest  and  attempts  were 
made  to  remove  the  mean  current  from  the  records  so  that  the 
magnitude,  damping  period,  and  associated  energy  of  the 
inertial  motion  could  be  determined.  In  the  first  attempt, 
the  mean  current  over  the  first  three  inertial  periods 
(pre-storra  passage)  was  subtracted  from  all  of  the  a  veloc- 
ity components  and  the  v  velocity  components.  The  three 
inertial  period  average  started  at  1920  GMT  07  September, 
which  is  five  days  before  the  storm  center  passed  the  near- 
est buoy.  Comparison  of  the  initial  PVD*s  (Fig.  11)  and 
similar  diagrams  of  the  same  lata  with  the  three  inertia! 
period  average  removed    (Fig.    12)     was    maie. 

The  initial  PVD's  indicate  a  general  transport  in  direc- 
tions  which,      in   most   cases,      are   nearly   the    same  before   and 
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Figure   1  1. 


Initial   PVD's    for    7-18   September    at    depths    of 
[a\     21    m    and    j[  b)     251     in  at   station    3.      Tha 

initial    location    is    at    (0,0). 
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Figure   1  2. 


PVD's  resulting  from    the    removal    of   the   ore- 
storm   averages   over    three   inertial   periods    from 
the    PVD'  s    of   Fig.    11  . 
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after  the  storm-induced  perturbation.  The  perturbations  are 
generally  reminiscent  of  inertial  motion,  i.e.  oscillations 
with  a  period  equal  to  the  local  inertial  period.  PVD's  of 
purely  inertial  motion  would  transcribe  a  simple  circle  for 
each  inertial  period.  The  circles  would  overlap  if  more 
than  one      inertial   period     were    plotted.  Where   the      first 

three  days  of  a  record  were  similar  to  the  mean  flow,  remov- 
ing the  three  inertial  period  average  from  the  record 
resulted  in  mere  nearly  overlapping  circular  PVD's.  Remov- 
ing the  three  inertial  period  average  was  not  effective  in 
isolating  the  inertial  component  of  the  flow  when  the  first 
three  days  of  the  record  were  not  obviously  in  the  general 
transport  direction.  Therefore,  the  mean  current  was  chang- 
ing significantly  during/following  the  storm  at  many  of  the 
current  meter  locations.  A  mora  effective  method  of  isolat- 
ing the  inertial  component  from  each  of  the  records  was 
required. 

D.        ISOLATING    THE    INERTIAL    MOTION 

In  a  further  attempt  to  isolate  the  large  inertial  com- 
ponent that  was  obvious  in  the  original  records,  a  running 
average  over  the  inertial  period  was  computed  each  three 
hours   and   then  subtracted    from  the   initial    ?VD   records.      The 
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averaging  technigue  is  necessarily  applied  at  the  central 
point  of  each  inertial  period,  and  as  a  result  the  half- 
periods  at  the  beginning  and  and  of  the  records  were  lost. 
However,  the  transient  response  is  not  great  at  these  times 
and   the      information   lost      is   inconsequential.  Next,      the 

average  inertial  period  record  is  linearly  interpolated 
between  3-h  values  to  obtain  the  same  nunber  of  data  points 
as  the  initial  P7D  data.  PVD's  of  the  running  inertial 
period  average  data  (e.g^.  Fig.  13)  show  that  the  majority  of 
the  inertial  flow  in  the  raw  data  is  removed  by  the  running 
average.  The  running  inertial  period  averaged  PVD's  show 
the  general  transport  direction  as  well  as  some  features 
which  should  be  noticeable  in  a  model  of  the  ocean  response. 
k  surge  to  the  northeast,  onshore,  is  noticeable  in  the 
record  of  Fig.  13a,  and  there  is  some  indication  that  the 
post-storm  current  is  greater  than  pre-storra,  i.e.  the  hori- 
zontal displacement  between  inertial  period  marks  is  greater 
after  the     storm   surge      than   before.  The    best      example   of 

these  features  in  any  of  the  records  is  shown  in  Fig.  13b. 
k  surge  to  the  west  followed  by  a  counterclockwise  loop  back 
into  the  pre-storm  flow  direction  is  evident.  Additionally, 
the    post-storm      flow   is   of      much    greater   magnitude     than    tne 
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pre-storm  flow.  These  diagrams  were  obviously  a  better  rep- 
resentation of  the  mean  flow  throughout  the  record  than  the 
previously  computed  three  inertial  period  averages.  These 
PVD's  also  show  that  the  mean  current  was  indeed  influenced 
by  the  bottom  topography.  Comparison  of  Fig.  4  and  repre- 
sentative PVD's  at  stations  one,  two  and  three  (Figs.  14  and 
15)  show  that  the  mean  flow  follows  the  bathymetric  con- 
tours. The  current  moves  rapidly  toward  the  northwest  a* 
station  one  during  and  for  some  time  following  the  period  of 
hurricane  passage  (Fig.  14a).  I  counterclockwise  loop  which 
would  be  associated  with  inertial  motion  is  not  seen  as  it 
is  in  the  station  two  record  (Fig.  14b)  .  Comparison  of 
Figs.  13b  and  15  shows  the  flow  is  along  the  contours  but  in 
opposite  directions  at  the  251  and  the  437  or  457  meter 
depths.  Also,  a  shoreward  surge  similar  to  that  at  the 
near-surface  level  of  station  three  (Fig.  13a)  is  seen  at 
least  one  inertial  period  later  near  the  bottom  (Fig.  15). 
The  tendency  to  return  to  the  pre-storm  flow  direction  is 
much  slower  at  the  near-bottom  records  and  is  not  oomple-e 
by    the  end   of  the   record. 

The   running    inertial-period    averages   are   then   subtracted 
from    the      initial   PVD    data      point-by-point    and   PVD's      of   the 
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Running   inertial   period   averagad    PyD   at   a    depth 
of    (a)    437    m   and    {b)     457    n    at    station    3. 
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resultant  data  are  plotted  (Fig.  16)  .  Examination  of  these 
resultant  PVD*s  shows  that  station  thrae  exhibits  the  most 
obvious  inertial  response,  possibly  dus  to  less  coastal  or 
bottom  influence.  The  resultant  PVD*s  for  stations  one  and 
two  (not  shown)  also  show  large  inertial  responses.  The 
station  one  record  still  exhibits  a  strong  flow  which  is 
non-inertial.  The      station      two      recori      shows     that      the 

majority    of  the    non-inertial  motion   is   removed. 

It  is  not  possible  from  these  plots  to  determine  pre- 
cisely when  the  response  to  the  hurricane  forcing  first 
begins  at  each  depth.  Since  the  rate  at  which  the  energy  is 
transmitted  vertically  is  of  interest,  tha  next  procedure  is 
to  calculate  the  energy  associated  with  each  inertial  period 
at    all  available   levels. 

E.        PROCEDURES    FOR    DETERMINING    RELATIVE    ENERGIES 

The  PVD's  which  were  computed  in  the  above  steps  indi- 
cate the  average  of  the  current  record  is  not  zero.  Also,  a 
breakdown  of  these  records  into  inertial  periods  indicate 
this  nonzero  mean  current  which  does  exist  is  not  constant 
throughout  the  record.  An  interactive  computer  program  was 
developed  to  remove  the  mean  current  from  a  record  consist- 
ing   of  inertial   flow    superposed    on   a   mean    current.      The   mean 
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velocity  is  calculated  as  the  distance  between  the  starting 
and  ending  points  of  the  record  divided  by  the  duration  of 
the  record.  An  egual  portion  of  this  mean  current  is  sub- 
tracted from  each  ten  minute  interval  of  the  record.  The 
mean  current  is  removed  from  the  record  one  inertial  period 
at  a  time  rather  than  for  an  entire  record  since  the  mean 
current  is  not  constant  throughout  the  entire  record. 

A  circle  of  radius  R  was  fit  to  each  of  the  inertial 
period  records.  The  circle  coordinates  (x,y)  are  computed 
using  a  mean  velocity.   Given  that 

X     =    R    COS  {0) 

y   =  R   sin(©> 
The    mean    velocity   is   computed   as      the   time    derivative    of   the 
position   as: 

u   =   dx/dt    =  -R   sin  «9)     d(9/dt 
v   =    dy/dt    =    R    zos(0)     d@/it 
where 

ae/at  =  -2  7t/(n  dt) 

and  N  represents  the  number  of  points  used  to  make  the  cir- 
cle. Successive  x  and  y  coordinates  of  the  circle  are  then 
computed    as: 

x1    =    xO    *     (u  dt) 

y1    =    yQ    +     (v   dt) 
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The  starting  position  is  x  ■  0.  and  y  =  -1.,  i.e.  at  270 
degrees  of  a  360  degree  circle  centered  at  (0.,0.),  and  the 
successive  positions  are  created  in  a  clockwise  fashion. 
The  radius  of  the  circle  is  adjusted  using  an  error  mini- 
mization program.  This  program  calculates  the  distances 
between  the  positions  reguired  to  create  the  circle  and  the 
positions  used  to  compute  the  PVD's,  and  selects  tha  circle 
which  gives  the  smallest  sum  of  positional  distances.  One 
fourth  of  the  maximum  dimension  of  the  plotted  data  is  usad 
as  the  initial  radius  for  the  minimization  subroutine. 
While  the  radius  is  held  constant,  the  sun  of  the  difference 
between  an  "average  error"  and  an  "individual  position 
error"  is  computed  and  divided  by  the  number  of  differences. 
This  error  difference  is  used  in  the  following  calculations. 
The  distances  between  all  points  of  the  computed  circle  and 
the  data  are  first  compared  to  locate  the  data  point  and  the 
circle  point  which  are  closest.  These  points  are  used  as 
the  starting  positions  for  the  distance  arror  computations. 
The  "average  error"  is  the  average  of  the  distances  between 
the  consecutive  points  of  the  data  and  the  computed  circle. 
The  "individual  position  error"  is  the  distance  between  a 
point  on  the  circle  and  a  data  point.    The  error  difference 
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for  one  radius  is  stored  and,  the  radius  is  slowly  incre- 
mented to  twice  the  original  magnitude,  storing  the  error 
difference  at  each  increment.  The  radius  corresponding  to 
the  minimum  value  of  error  difference  is  used  to  create  the 
circle  which  best  fits  the  plotted  data. 

The  circle  created  using  this  adjusted  radius  is  used  as 
the  best  fit.  It  is  then  possible  to  estimate  the  energy 
per  unit  mass  associated  with  each  inertial  period.  The 
inertial  velocity  (V=fH)  associated  with  the  circle  radius  H 
and  coriolis  parameter  f  is  computed.  One  half  of  the 
square  of  this  velocity  is  assumed  to  be  the  energy  per  unit 
mass.  The  circle  fitting  error  difference  computed  as  above 
is  divided  by  the  radius  of  the  fitted  circle  and  referred 
to  as  the  relative  RMS  error.  The  energy  values  are  plotted 
versus  the  inertial  period  at  each  level  at  each  station 
(See  Figs.  17  and  18).  The  radius  for  each  inertial  period, 
the  relative  RMS  error  of  fitting  the  circle  and  the  mean 
current  removed  from  each  inertial  period  are  listed  in 
Table  II  for  the  records  from  station  one  at  a  depth  of  49  m 
and  from  station  three  at  a  depth  of  21  m.  Also,  the  mini- 
mization program  is  tested  using  exact  circles  as  input 
data.   The  resultant  errors  for  given  radii  are  also  listed 
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in  Table  II.  The  exact  circle  error  shows  the  error  is  near 
zero  when  the  motion  is  circular.  The  error  is  found  to  be 
a  smaller  percentage  of  the  radius,  i.  e.  the  relative  RMS 
error  is  smaller,  at  station  three  than  at  station  one. 
This  supports  the  previous  results  of  the  motion  being  more 
nearly  circular,  and  thus  mora  inertial,  at  station  three 
than  at  station  one. 

A  good  fit  between  the  positions  of  tie  generated  circle 
and  the  positions  from  the  real  data  with  the  inertial  mean 
removed  is  obtained  for  the  records  of  station  three  where 
the  water  is  deepest  and  the  coastal  influence  is  least 
(Fig.  19a).  The  agreement  in  positions  between  the  two 
curves  at  station  one  is  generally  less  since  these  records 
are  affected  more  by  the  bottom,  and  by  tidal  motion  which 
is  not  necessarily  removed.  Obviously  non-inertial  motion 
remained  in  the  station  one  record  (Fig.  19b).  Example 
plots  of  the  record  after  all  means  were  removed  are  shown 
in  Fig.  20.  The  figures  show  nearly  overlapping  circular 
motion  except  for  different  radii.  Plots  of  the  energy  val- 
ues versus  inertial  periods  are  given  in  Fig.  17  for  the 
level  nearest  the  surface  at  each  of  the  three  stations. 
Energy   values   of   4000   at  station   3   correspond   to   an 
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TABLE  II 

Example  circle  radii,  relative  RMS  error  and  mean  current 
removed  for  each  mertial  period  ana  for  an  exact  circle 


RELATIVE 

INERTIAL 

RMS 

CIRCLE 

CURRE 

STATION 

PERIOD 

ERROR 

RADIUS 

REMOVED 

(km) 

(Ha) 

(cm 
u~' 

sii) 

V 

1 

1 

.3505 

.321 

-    .  13 

.04 

2 

.3775 

.631 

-    .39 

-    .56 

3 

.3810 

.839 

-5.18 

.38 

4 

.2622 

1.  152 

1.15 

1.69 

5 

.1663 

2.663 

-6.53 

1.33 

6 

.1688 

2.602 

-    .52 

2.03 

7 

.1370 

2.859 

6.18 

-10.77 

8 

.1580 

2.577 

1.22 

2.44 

9 

.2110 

1.123 

.24 

.43 

10 

.2219 

1.560 

-3.23 

4.73 

11 

.2165 

1.714 

1.  16 

-2.62 

12 

.2695 

1.302 

2.89 

-1.64 

3 

1 

.0818 

2.547 

-    .01 

.29 

2 

.0738 

2.335 

1.83 

-    .  29 

3 

.1501 

2.106 

.19 

-1.  17 

4 

.1530 

2.035 

.20 

-    .92 

5 

.2813 

1.440 

-5.67 

2.92 

6 

.0339 

12.455 

-9.66 

-2.81 

7 

.0322 

12.116 

-1.83 

5.89 

8 

.0160 

13.394 

3.52 

4.75 

9 

.0316 

1 0.904 

3.56 

-    .16 

10 

.0301 

6.357 

2.07 

3.77 

11 

.0583 

5.818 

-    .67 

.  89 

12 

.0853 

3.891 

-2.00 

-2.  16 

True   Ci 

.rcle 

.0033 
.0022 
.0002 
.0003 
.0002 
.0005 
.0023 
.0002 
.0021 

.300 

.500 

1.000 

2.000 

3.000 

4.000 

8.000 

16.000 

3  0.000 

amplitude  of  the  inertial  curreat  of  about  90  cm  s~l.  The 
energy  values  at  station  one  are  much  smaller  than  those  at 
station  two  and  three.    The  technique  employed   to  fit  the 
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true  circle  to  the  manipulated  data  accounts  for  some  of 
this  difference.  When  the  data  points  are  not  circular,  the 
circle  is  adjusted  to  some  value  less  than  the  largest 
dimension  of  the  plotted  data.  When  the  data  form  a  circle, 
the  largest  dimension  of  the  plotted  data  and  the  dimension 
of  the  circle  nearly  coincide.  For  stations  two  and  three, 
a  sharp  increase  in  kinetic  energy  is  indicated  in  Fig.  17 
during      the   fifth      inertial     period.  These   energy      values 

remain  high  for  two  to  three  iner+ial  periods  before  rapidly 
decreasing  toward  the  pre-storm  levels.  Similar  plots  for 
the  next  deepest  level  at  each  station  are  shown  in 
Fig.  18a.  At  stations  one  and  two  the  energy  peaks  rapidly, 
falls  off  rapidly  and  then  increases  to  the  end  of  the 
record.  At  station  three,  a  slower  increase  is  noticed,  and 
is  followed  by  a  nearly  constant  value  to  the  end  of  the 
record.  A  slower  increase  in  energy  is  found  (Fig.  18b)  at 
the    437      n   and   457   m      levels   at    station   three.  The    energy 

values  do  not  appear  to  have  peaked  by  the  end  of  the 
record.  Since  the  records  of  Fig.  18  are  within  20  m  of 
each  other,  it  would  be  reasonable  to  assume  the  records 
would  be  nearly  identical.  The  difference  could  be  due  to 
bottom  reflection   of      energy  since   the   '457      m   instrument   was 
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positioned  only  8  meters  from  the  bottom.  At  the  very 
least,  the  procedures  used  to  determine  the  energy  associ- 
ated with  the  inertial  motion  should  produce  results  accu- 
rate to  within  a  factor  of  two  and  should  also  indicate  the 
trends.  No  confidence  is  place!  in  the  details  of  tae  fluc- 
tuating values.  However,  the  large  peaks  and  long  term 
increases  shown  in  the  energy  records  are  thought  to  be 
significant. 
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Figure  17. 


Energy  (m2  s~2)  versus  inertial  period  at  the 
near  surface  levels  at  stations  one,  two  and 
three. 
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Energy, (m2   sr2)    versus   lnertial    Period   at     (a) 
the    mid    levels   of   aach  station    and    (b)    the    two 
bottom    levels    at   station    threa. 
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Manipulated   real   data    record   and    calculated 
circle    for    (a),    inert  lal   penol    6    at   station   3 
and     (b)    inertial   period   6    at    station    1. 
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III.       THE    EMBEDDED    MIXED   LAYER--OCEAN 
CTFCUUTlOff-TfODEI 


A.        EACK GROUND 

Adamec  et  al  (1981)  tested  an  embsdded  mixed  layer- 
ocean  circulation  model.  This  model  ambedded  Garwood's 
(1977)  model  for  predicting  mixed  layer  depth  and  jumps  in 
temperature  and  velocity  at  the  base  of  the  mixed  layer  in  a 
six-level,  primitive  equation  open-ocean  model  developed  by 
Haney  (198  1).  These  authors  remarked  in  their  conclusion 
that  additional  tests,  "including  comparisons  with  observa- 
tional data",  were  being  planned.  This  thesis  constitutes 
an  important  step  in  those  plans  because  of  the  availability 
of  suitable  observational  data.  The  model  used  in  the  1931 
study  is  expanded  to  simulate  the  thres-dimensional  ocean 
response  to  a  translating  hurricane-type  forcing  field. 
This  model  forcing  represents  an  axisymmetric  hurricane 
translating  at  a  constant  speed  of  7.5  a  s~ l .  The  model 
ocean  basin  is  a  960  km  square  with  six  layers  of  varying 
depth  between  the  surface  and  a  free-slip  bottom  at  400  m. 
The  horizontal  resolution  is  15  km  on  a  65  by  65  grid.  The 
level  depths  and  thicknesses  are  represented  in  Fig.  21. 
The    time    step   used    in   the    198  1    study,    450    s,    is    increased    to 
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600    s.  A   brief  explanation   of      the   modal   as   used      in   this 

thesis  is   given    below. 

B.        MODEL    FORMULATIONS    AND    BOUNDARY    CONDITIONS 

The  governing  equations  ara  the  primitive  Navier-Stokes 
equations,  the  continuity  equation,  the  hydrostatic  equation 
and  the  equation  of  state  (sea  Adamec  at  al,  1981).  The 
hydrcstatic  assumption  is  appliad,  and  the  ocean  is  assumed 
incompressible  with  density  being  a  linear  function  of  temp- 
erature only.  The  coriolis  force  varias  appropriate  to  the 
latitude  of  the  domain  from  25JN  to  33.6°N.  The  change  in 
coriolis  parameter  (f)  with  latituda  is  computed  as  a  finite 
difference    at   each  north-south  grid   point. 

There  are  no  fluxes  of  mass,  momentum  or  heat  across  the 
bottom  (flat)  or  side  (vertical)  boundaries.  Also,  the 
rigid  lid  approximation  (w  ■  0  at  z  ■  0)  is  made,  which 
requires  the  vertically  averaged  motion  be  zero.  Applying 
this  approximation  to  the  Navia r-Stokes  equations  provides 
prediction  equations  for  the  vertical  shear  currents  (see 
Adamec  et    al,       1981).  Vertically   averaging    the   continuity 

equation  over  the  mixed  layer  and  applying  w  =  0  at  z  -  0 
yields  a  prognostic  aquation  for  the  vertical  velocity  at 
the    base    of   the    well-mixed      layer,      w  (-h) .         The    entrainment 
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velocity  is  predicted  by  the  mixed  layer  model.  Together 
with  w(-h),  this  yields  a  prognostic  aquation  for  mixed 
layer  depth. 

The  axisymmetrical  portions  of  the  hurricane  forcing  are 
shown  in  Fig.  22.  The  tangential  and  radial  stresses  are  of 
the  same  form.  They  are  calculated  from  a  prescribed  wind, 
which  increases  linearly  from  the  storm  eye  wall  to  the 
radius  of  maximum  winds,  and  then  decreases  as  r~l/2  to  r  = 
360  km  and  then  linearly  to  r  =  450  Jem.  The  radius  of  maxi- 
mum winds  is  45  km  and  the  inner  boundary  of  the  eye  wall  of 
the  storm  is  4.5  km.  This  wind  profile  results  in  a  wind 
stress  curl  which  is  zero  inside  the  eye  wall,  increasing 
linearly  from  zero  at  the  eye  wall  to  a  maximum  at  45  km, 
zero  from  45  to  360  km,  and  negative  from  360  km  to  the 
boundary  of  the  wind  stress  (about  450  km)  .  The  maximum 
tangential  stress  corresponds  to  a  maximum  wind  speed  of  50 
m  s-1.  The  small  value  of  radial  stress  (-12.9  iPa  max)  is 
due  to  cross  isobaric  flow  of  about  20°. 

Eelow  the  mixed  layer,  vertical  diffusion  is  applied  to 
the  momentum  or  temperature  equations  with  a  vertical  eddy 
viscosity  or  eddy  conductivity  coefficenr  (both  equal  to  0.5 
cm  s-2  at  all  depths  below  the  surface  mixed  layer).   In  the 
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Radius  [km] 

Figure  22.      Forcing   functions   for   the   ocean    model,    each 

normalized   by   the   value   of   the    function   at    the 
radius   of   maximum  winds.      The    tangential    and 
radial   stress    components   have   the   same   form 
(solid  line)    with   maximum   values   of   35.9   and 
-12.9   dPa.    respectively.      The    surface   heat    flux 
(dashed)     has   a    maximum  value    of   -840    W   m 
(Adamec  et   al,    1981)  . 

mixed     layer,         however,        the        above      formulation      is      not 

appropriate.        The  depth   of  the      mixed   layer   is   the  boundary 

between      intense      turbulence     and   the      much      less     turbulent 

waters  beneath. 


C.        ENTRAINMENT    AND    MIXED    LAYER     MODEL    FORMULATION 

In  the  mixed  layer,  prognostic  equations  for  the  mixed- 
layer  average  (bulk)  values  of  the  vertical  component  of 
turbulent  kinetic  energy  and  the  total  turbulent  kinetic 
energy  are  derived  using  the  bulk,  seoond-order  closure 
methods  of   Garwood    (1977).  Computation    of    the   entrainment 

buoyancy    flux   allows      calculation    of    the    iownward      fluxes    of 
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heat  and  momentum  associated  with  entrainmer.t  at  the  base  of 
the  mixed  layer. 

The  above  mixed- layer  formulation  allows  the  fraction  of 
turbulent  kinetic  energy  available  for  mixing  at  the  base  of 
the  mixed  layer  to  be  other  than  a  constant.  This  turbulent 
kinetic  energy  is  diagnost ically  dependent  upon  the  surface 
buoyancy  flux  and  the  surface  friction  velocity  (known 
boundary  conditions) .  It  is  thus  possible  to  couple  the 
diagnostic  mixed-layer  formulation  with  the  prognostic  ocean 
circulation  model  in  a  numerically  feasible  fashion. 

Eoth  increases  and  decreases  in  mixed  layer  depth  must 
be  considered  in  any  mixed  layer  model.  The  easier  of  these 
two  events  to  formulate  is  the  increasing  mixed  layer  depth 
case.  The  method  requires  determination  of  the  entrainment 
heat  flux  and  then  imposing  this  heat  flux  on  the  given 
temperature  profile.  Added  vertical  resolution  near  the 
base  of  the  mixed  layer  is  given  since  the  base  is  not 
required  to  coincide  with  any  of  the  prescribed  model  lev- 
els. This  increased  vertical  resolution  is  very  important 
since  the  thermocline  profile  determines  the  potential 
energy  of  the  upper  ocean.  The  dynamics  of  the  mixed  layer 
and  the  ocean  circulation  are  dependent  on  this  potential 
energy. 
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Formation  of  a  new,  shallower  mixed  layer  occurs 
henever  the  turbulent  mixing  is  not  abla  to  penetrate  all 
the  way  to  the  previously  established  depth  of  the  mixed 
layer.  This  layer  reformation  occurs  when  warming  of  the 
surface  layer  occurs  with  no  increase  in  the  wind  shear  to 
sustain  the  mixing.  Numerical  formulation  of  this  event  is 
difficult  because  the  previous  structure  at  the  interface 
can  not  readily  be  preserved.  Thus  energy  budget  problems 
may  arise  when  deepening  to  the  prior  mixed  layer  depth  does 
occur.  This  model  uses  a  numerical  procedure  (see  Adamec  et 
al,  198  1)  that  preserves  potential  energy  to  ensure  the 
deepening  rate  is  as  correct  as  possible  when  the  layer 
again  deepens  to  the  earlier  mixed  layer  depth.  This  fea- 
ture of  the  model  is  not  tested  here  because  the  solar  flux 
is  zero  throughout  the  simulation. 

D.   DYNAMIC  STABILITY  CONDITION 

It  is  assumed  that  the  mixed  layer  is  dynamically  unsta- 
ble and  the  underlying  water  column  is  normally  dynamically 
stable.  However,  dynamic  instability  of  the  underlying 
water  column  can  sometimes  occur.  In  this  model,  vertical 
fluxes  of  heat  and  momentum  between  levels  are  imposed  so 
that  the  gradient  Richardson  number   remains  greater  than  cr 
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equal  to  a  critical  value-  This  generalization  of  the  more 
common  convective  adjustment  is  referred  to  by  Adamec  et  al 
(198  1)    as    "dynamical   adjustment". 

E.  COUPLING    OF    THE    DYNAMICAL    AND    MIXING    PROCESSES    IN    THE 
MODEL 

The  mixed  layer  and  the  dynamic  portion  of  the  model, 
which  is  a  level  model  that  predicts  the  average  of  a  quan- 
tity in  a  layer,  are  coupled  in  two  phases.  First,  advec- 
tive  and  diffusive  changes  in  the  upper  ocean  are  calculated 
in  the  dynamic  part  and  put  into  a  form  useable  by  the  mixed 
layer  model.  Then,  the  surface  flux  and  sntrainment  changes 
are  calculated  by  the  mixed  layer  model  and  transmitted  to 
the    dynamic   part   of   the   model.  A   special    treatment    of   the 

level   which  contains   the   base   of    the    mixed    layer   is  required 
(see   Adamec  et   al,    1981)  . 

F.  COMPARISON  OF  THE  MODEL  AND  THE  DATA 

The  model  provides  a  four-dimensional  space-time  view  of 
the  ocean  response  to  a  moving  hurricane.  The  data  provided 
by  NAVOCEANO  is  two-dimensional,  time  and  depth,  at  three 
different  stations.  All  three  locations  are  to  the  right  of 
the  storm  track  beyond  the  radius  of  maximum  winds  of  the 
storm.      The   differences    in    instrument   depths   and   the    limited 
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number  of  stations  do  not  allow  the  extension  of  the  data 
into  a  third  dimension.  The  comparison  of  the  data  and  the 
model,   then,    is   an  attempt    to: 

•  compare  the  general  physical  characteristics  of  the 
model  results  to  see  if  they  are  reasonable  in  light  of 
the    observed   data   results, 

•  choose  a  data  set  from  the  model  which  is  similar  to 
the  raw  data  in  number,  depth  and  location  relative  to 
the   storm, 

•  compare  this  data  to  the  NAVOCEANO  data  to  show  whether 
or  not  the  model  is  providing  a  realistic  response  to 
the    simulated   hurricane   forcing,    and 

•  examine  the  four-dimensional  model  results  for  charac- 
teristic ocean  response  to  hurricane  forcing  which 
could  not   be   observed   at   the  data   stations. 

G.        GENERAL    CHARACTERISTICS   OF    THE    MODEL 

The   initial    position      of   the    hurricane    is      at   grid    point 
(23,16).  Since   the      grid  spacing      is   15      km,      grid      point 

(23,16)  corresponds  to  a  distance  in  the  x  direction  of  33C 
km  and  in  the  y  direction  of  225  km.  Moiel  output  includes 
north-south  and  east-west  velocity  components,  v  and  u  ; 
temperature,  T  ;  and  mixed  layer  depth,  H.  The  initial  val- 
ues cf  these  guanties  are  given  in  the  following  table. 
The  guantities  u,  v,  and  T  are  extracted  at  both  three  hour 
intervals  for  the  entire  grid  at  all  six  levels  and  at  ten 
minute  intervals  for  selected  grid  points  and  levels.  The 
mixed   layer  depth    (MLD)       is     stored  at   three-hour   intervals. 
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TABLE   III 
Initial   values   of  model   output  variables 

I    at   level 
u        Z  HI  2  3  4  5  6 

(cm   s-i)  (m)  (°C) 

0        0  30.         30.0        29.2         28.1         23.4         18.4         13.1 

The  results  predicted  by  the  model  at  36  and  43  h  at  level  1 
are  presented  in  Figs.  23a--d  and  Figs.  24a — d,  respec- 
tively. Using  the  storm  translation  speed  of  7.5  m  s~lf  the 
distance  of  the  storm  from  the  bottom  of  the  grid  is  given 
in  Table  IV.  The  storm  center  has  tracked  off  the  grid  by 
hour  24  and  the  southern  boundary  of  the  storm  has  passed 
the    edge    by  hour   48. 

TABLE    IV 
Hurricane    distance   from    bottom   of  grid   at    six-hour   intervals 

DISTANCE  HOUR 

330  0 

492  6 

654  12 

816  18 

978  24 

1140  30 

1302  36 

1464  42 

Several   important      features    are     noticeable   in     Figs.    23 

and    24: 

•  The  current  fields  at  36  and  43  hours  are  "out  of 
phase",  i.e.  the  areas  of,  westerly  currents  at  hour  36 
are    generally    easterly  at   hour  48. 
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Model   output    of    (a)     u   velocity    component,     (b)     ' 
velocity    comoonent,     (c)    mixsd-layer    deptn,    and 
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•  The  resultant  motion  is  circular  over  nearly  the  entire 
field  with  an  apparent  rotation  period  of  about  24 
hours. 

•  Current  velocities  in  the  mixed  layer  range  from  neg- 
ative values  of  about  70  cm  - l  to  positive  values  of 
about   80  cm   s~l. 

•  The  mixed  layer  deepens  preferentially  to  the  right  of 
the  storm,  and  ranges  from  minimum  values  of  about  30  m 
to   maximum    values   of   almost    60  m. 

•  A  ridge  of  cold,  upwelled  water  lies  in  the  wake  of  the 
storm.  This  ridge  is  shown  best  in  the  level  3  temper- 
ature contour  of  Fig.  25.  Superposed  on  this  cold 
ridge  is  an  upwelling/downwelling  wake  with  an  inertial 
period.  Upwelling  is  in  the  areas  of  minimum  MLD's 
(32.5  in  Fig.  24c) ,  and  downwelling  is  in  the  area  of 
maximum  MLD's  (58.0  in  Fig.  23c).  Minimum  MLD's  are 
reached  at  the  end  of  the  upwelling  cycle  and  maximum 
MLD's   are  reached   at   the    end  of   the   iownwelling   cycle. 

These  results  are  consistent  with  the  results  derived  from 
the  NAVOCEANO  data,  which  shows  that  the  primary  ocean 
current  response  to  hurricane  passage  when  no  mean  current 
exists  is  circular  motion  at  the  local  inertial  period.  The 
reversal  in  flow  direction  every  12  hours  is  particularly 
evident  in  the  u  velocity  component  records  for  level  six 
shown  in  Fig.  26.  The  model  predictions  describe  the  ocean 
response  over  the  entire  grid,  rather  than  just  at  some 
selected  points.  Trajectories  are  plotted  at  several  loca- 
tions perpendicular  to  the  storm  track  (see  Fig.  27) .  The 
surface-layer  trajectory  that  began  on  the  left  side  of  the 
track  (Figs.  27a)  clearly  shows  the  surge  to  the  left  as 
the    storm   approaches,    strong   inertial   motion    following   storm 


65 


passage  and  a  net  transport  to  the  south.  The  trajectory 
along  the  storm  path  (Fig.  27b)  shows  a  surge  to  the  left, 
recovery  toward  the  right  and  a  small  net  displacement.  The 
trajectories  that  began  on  the  right  (Figs.  27c  and  27d) 
show  a  net  deflection  northward  and  to  the  right.  These 
trajectories  show  that  different  transport  profiles  are 
expected  at  different  locations  relative  to  the  stora.  This 
feature  is  not  obvious  in  the  NAVOCEAND  data  due  to  the 
limited  number  of  locations  of  the  observations. 

B.   SELECTED  POSITION  AND  DEPTH  DATA 

A  data  set  similar  to  the  NAVOCEANO  data  set  was  derived 
by  extracting  model  variables  at  ten  minute  intervals  at  the 
three  positions  and  depths  shown  in  Table  V.  These  time 
series  of  model  variables  (Figs.  28  and  29)  are  plotted  sim- 
ilarly to  the  raw  data  plots  of  Figs.  6  and  7.  The  observed 
(Figs.  6)  and  predicted  (Fig.  28)  u-::o:nponents  show  very 
similar  variations  although  the  model  data  are  much 
smoother.  Also,  although  the  magnitude  of  the  velocity  com- 
ponent at  level  2  of  Fig.  28a  is  somewhat  smaller  than  the 
21  m  magnitude  shown  in  Fig.  5a,  the  magnitude  at  level  6 
(300  m)  of  Fig.  28b  is  much  smaller  than  the  437  m  magni- 
tude shown  in   Fig.  6b.    The  implication  is  that  not  enough 
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Model   temperature   prediction   for    level   three   a* 
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Figure   27. 


Trajectories   for    level   1    at   Doints    (a)    60    km 
left   of,     (b)    along,     (c)    60    km    right   of,    (d)     120 
km   right   of    storm   track.       Each    tick   mark   along 
the   trajectory    represents   a   3-h   displacement.      X 
indicates   the    initial   position    of   each 
trajectory.      Notice    each    plot   employs   different 
spatial   increments. 
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energy  is  being  vertically  transmitted  by  the  model.  Both 
the  model  motion  and  the  observed  flow  are  inertial  and  nei- 
ther has  completely  damped  by  the  end  of  nine  days.  The 
model  a  components  of  the  velocity  at  levels  one  and  six  are 
exactly  180  degrees  out  of  phase  (Fig.  30a) .  That  is,  the 
velocities  of  the  surface  and  bottom  model  layers  are  in  the 
opposite  directions.  The  NAVOCEANO  velocities  are  also 
exactly  out  of  phase  immediately  after  storm  passage. 
although  they  become  more  in  phase  later  in  the  record  (see 
Fig.  30b) .  For  any  given  level  of  the  nodel,  the  u  and  v 
variations   are  90    degrees   cut   of    phase. 

TABLE    V 

Positions    and  depths    for   model   ten  minute   interval   data   and 

corresponding   NAVOCEANO   data. 

DISTANCE 

INERTIAL                           FROM  MID  NAVOCEANO 

GRID                      PERIOD                    STORM    CENTER  LAYER  INSTRUMENT 

POSITION    MODEL    NAVOCEANO       MODEL     NAVOCEANO  DEPTH  DEPTH 

(x,y)                   (hours)                          (km)  (m)  (m) 

30,33  2U.50         24.54 

29,35  24.30         24.34 

31,39  23.95         23.90 

Three-hour   running   averages    by  inertial    period   were   com- 
puted  for    the   extracted    model      data   using    the    same  technique 
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Superpositions   of   the  surface   and   deep   layer 
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as  used  for  the  NAVOCEANO  results.  A  direct  comparison  of 
Figs.  13,  14  and  15  to  the  extracted  model  results  running 
inertial  period  averages  is  mada  in  Figs.  31,  32,  33.  The 
model  near-surface  residual  flows  are  very  similar  to  each 
other  (Figs.  31a  and  32a)  ,  as  are  the  lower  level  results 
(Figs.  31b,  32b  and  33).  Net  transport  to  the  north  is 
indicated  in  all  records  of  the  modal.  The  NAVOCEANO  data 
indicate  net  transports  that  are  closely  tied  to  the  local 
topography,      which   is      absent   in    the   model.  The  NAVOCEANO 

results  also  show  opposite  current  directions  for  near-sur- 
face and  mid  depth.  These  observations  provide  interesting 
examples  of  results  not  shown  in  the  model,  and  which 
require   addtional   investigation    to  explain. 

In  order  to  make  a  more  quantitative  comparison  between 
the  model  and  the  NAVOCEANO  inertial  oscillations,  the  model 
results  from  the  selected  stations  shown  in  Table  V  were 
treated  in  the  same  manner  as  the  NAVOCEANO  raw  data.  Cir- 
cle radii,  relative  RMS  error  and  mean  current  removed  are 
listed    (Table      VII)      similarly   to     Table   II.  Relative   RMS 

errors  are  similar  to  those  for  station  three  (Table  II) 
although  maximum  circle  radii  are  not  as  large.  Resultant 
calculations   of    energy    (cm2  s-2)       versus   inertial   period   are 
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Figure  33. 


Running   inertial   period    average   PVD's   for    1( 
days   at    (a)     level   5    and    (b)     lavel   6    at   grid 
point    31,39. 
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plotted  (Figs.  34  and  35)  similarly  to  the  energy  plots  of 
Figs.  17  and  18.  These  figures  show  that  the  major  differ- 
ences in  the  model  and  the  NAVOCEANO  data,  at  least  in  a 
qualitative  sense  are  in  the  rate  of  propagation  of  energy 
with  depth.  It  appears  that  less  energy  is  propagated  ver- 
tically in  the  model  compared  to  the  observations.  The 
energies  in  Figs.  35a  and  35b  are  much  smaller  in  comparison 
to  the  values  in  Fig.  34  than  are  the  similar  results  of  the 
observed  data.  Also,  there  is  no  slow  increase  in  the  hori- 
zontal kinetic  energy  near  the  bottom  in  the  model  as 
observed  in  the  two  deepest  NA70CEAN0  records.  The  level 
four  records  in  Fig.  34a  suggest  that  energy  is  being  con- 
tinually propagated  into  the  layer  because  the  energy  levels 
do  net  appreciably  decay.  Detailed  calculations  are  neces- 
sary to  establish  the  energy  flux  mechanisms  from  the  model 
data. 


TABLE    VI 

Sample  of   circle   radii,    relative    RMS    error    and   mean   current 
removed   for   each   inertial    period    from   the   model   simulation. 


LEVEL 

RELATIVE 

ONE 

INERTIAL 

RMS 

CIRCLE 

CORRE 

STATION 

PERIOD 

ERROR 

RADIUS 

REMOVED 

(*,1) 

(*I) 

(kffi) 

(cm 
u 

§zi) 

V 

(30,33) 

1 

.0530 

6.373 

2.33 

2.73 

2 

.0334 

10.166 

-5.28 

5.63 
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IV.   SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 

A.   SUMMARY 

The  NAVOCEANO  data  provided  several  physically 
significant  results.  The  currents  measured  by  the 
instruments  were  affected  by  topography  and  by  the  storm. 
The  flow  tended  to  be  along  the  bathymetry  contours  except 
that  the  storm  imparted  a  surge  to  the  flow  consistent  with 
wind    stress     associated    with     the  storm.  The   station      one 

records  do  not  depict  easily  identifiable  inertial  motion. 
Station  one  is  in  shallow  water  of  about  100  meters,  and 
shoreward  motion  to  mid-depth  is  evident  prior  to  the  storm 
passage.  An  offshore  flow  is  evident  from  about  the  time  of 
the  storm  passage  for  a  duration  of  about  two  inertial  peri- 
ods. The  tendency  at  each  station  is  for  the  post-storm 
flow  to  return  to  the  same  direction  as  the  pre-s-orm  flow. 
Energy  input  by  the  storm  has  generally  dissipated  by  a  fac- 
tor of  e-1  after  three  or  four  inertial  periods,  and  the 
post-storm  flow  magnitude  is  greater  than  the  pre-storm  mag- 
nitude. The  horizontal  kinetic  energy  associated  with  iner- 
tial motion  in  the  near-surface  water  increased  more  than 
two  orders  of  magnitude  as  the  storm  passed  and  then  decayed 
over      six      or      seven      inertial      periods.  At      the      deepest 
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instrument  levels,  the  kinetic  energy  associated  with  iner- 
tial  motion  increased  slowly  and  had  not  peaked  seven  iner- 
tial  periods  after  storm  passage.  The  mid-level  instruments 
recorded  either  a  sharp  peak  after  storm  passage  or  a  less 
rapid  increase  and  little  decay  to  the  end  of  the  record 
depending  upon  location.  The  near-surface  and  near-bottom 
flows  are  about  180°  out  of  phase  for  at  least  five  or  six 
inertial    periods  after   storm   passage. 

The  mixed  layer-ocean  circulation  model  results  are  very 
similar  in  many  respects  to  the  NAVOCEANO  observations.  The 
model  does  not  include  bathymetry,  and  therefore  bathymetric 
effects  are  not  reproduced.  The  model  does  show,  however,  a 
storm-induced  surge,  inertial  motion  damping  over  several 
inertial  periods  at  the  surface  and  bottom  layers  and  dif- 
ferent magnitudes  in  pre-storm  and  post-storm  flow.  The 
model  predictions  indicate  an  almost  instantaneous  near-in- 
ertial  energy  propagation  rate  with  depth.  The  flow  at  the 
surface  and  at  depth  have  the  same  form,  increasing  in  mag- 
nitude and  decaying  at  very  nearly  the  same  rate  at  each 
level.  At  mid  levels  of  the  model,  the  kinetic  energy  level 
remains  nearly  constant.  The  current  oscillations  in  the 
mixed  layer    and      in    the   bottom   layer    are    180°      out   of    phase. 
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Trajectories  of  the  model  results  show  that  the  net  dis- 
placement at  some  stations  is  nearly  zero  while  at  other 
stations  the  net  displacement  is  tens  of  kilometers.  Also, 
the  direction  of  surge  and  of  the  residual  flow  is  dependent 
on  location  relative  to  the  storm  track.  The  model  shows  a 
zone  of  upwelling/downwelling  in  the  wake  of  the  storm  which 
has    a   naar-inertial   frequency. 

B.        CONCLUSIONS 

The  bathymetry  of  the  ocean  is  very  important  in  determ- 
ing  the  direction  of  flow  of  currents  at  all  levels.  This 
part  of  the  flow  is  not  included  in  this  ocean  circulation 
model  and  must  be  subtracted  from  the  real  ocean  data  if 
consistent  comparisons  of  model  results  and  current  observa- 
tions are  to  be  made.  The  effect  of  the  storm  is  to  drive 
the  average  motion  in  a  direction  consistent  with  the  wind 
field  of  the  storm.  Changes  in  direction  occur  over  a  very 
short  time  period  and  are  observed  at  all  depths.  As  the 
storm  passes,  inertial  motion  results  due  to  a  balance 
between   the      centrifugal  and     coriolis   accelerations.  The 

kinetic  energy  associated  with  this  motion  remains  nearly 
constant  near  the  surface  for  two  to  three  inertial  periods 
and    decays   over    an   additional   three   to   four    inertial   periods 
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to  a  post-storm  level  greater  than  the  pra-storm  level.  The 
available  data  do  not  show  the  energy  levels  returning  to 
the  pre- storm  level  before  the  end  of  the  record.  There  is 
an  obvious  time  lag  of  the  storm  which  increases  with  depth. 
The  exact  time  of  the  initial  affect  at  any  particular  depth 
is  difficult  to  measure  since  the  magnitude  of  the  velocity 
does  not  instantaneously  increase,  but  increases  slowly  with 
the  onset  of  the  storm.  Nonetheless,  tha  downward  propaga- 
tion of  this  energy  is  net  instantaneous  but  requires  on  the 
order  of  one-half  inertial  period  to  propagate  one-half  kil- 
ometer. The  energy  propagates  to  a  depth  of  about  450  m  in 
one-half  of  an  inertial  period.  At  this  depth  tha  energy 
increases  for  at  least  seven  to  eight  inertial  periods.  The 
mixed  layer-ocean  circulation  model  closely  simulates  the 
inertial  motion  and  damping  rates  near  the  surface.  The 
vertical  propagation  of  energy  is  not,  however,  very  similar 
to      the   NAVOCEANO      observations.  The      model    does      provide 

insight  into  the  three-dimensiD nal  ocean  and  is  impressive 
in  distinguishing  between  results  at  locations  varying  with 
respect   to    the   storm   center. 
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C.   BECOMMENDATIONS 

The  primary  differences  between  the  NAVOCEANO  data  and 
the  model  results  are  in  the  vertical  propagation  of  kinetic 
energy  associated  with  storm  passage.  Further  research  is 
required  to  isolate  the  mechanism  of  energy  propagation  so 
that  it  can  be  included  in  the  model  formulation. 
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